Periodic orbit theory in fractal drums 
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Abstract 

The level statistics of pseudointegrable fractal drums is studied numerically using 
periodic orbit theory. We find that the spectral rigidity As(L), which is a measure 
for the correlations between the eigenvalues, decreases to quite small values (as 
compared to systems with only small boundary roughness), thereby approaching 
the behavior of chaotic systems. The periodic orbit results are in good agreement 
with direct calculations of A%(L) from the eigenvalues. 
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1 Introduction 



Many systems with irregular geometries exist in nature and their physical 
properties (e.g. their electronic and vibrational behavior) have been widely 
investigated in the past. One class of irregular systems are those, which are 
formed by an ordered material, but possess an irregular shape of the boundary. 
A prominent example are the fractal drums [1], which are constructed by 
applying the so-called fractal generator (see Fig. 1(a)) v times to a regular 
square or rectangle. Already for prefractal shapes as e.g. for v — 1 — 3, the 
eigenstates possess energy spectra and localization properties that are well 
distinct from systems with smooth boundaries. This has been demonstrated by 
numerical simulations [2,3,4] as well as by experiments on liquid crystal films 
[4] and acoustic cavities [5]. In this paper, we address the question how the 
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specific geometry of the systems influence the behavior of the energy spectra 
within periodic-orbit theory, which provides a link between the classical and 
the quantum mechanical behavior of systems with the same shape. 

Two universality classes exist with a different classical dynamics, the chaotic 
and the regular systems. Fractal drums belong to the intermediate class of 
pseudointegrable systems, which are polygons with only rational angles nj7r/mj, 
with rii,mi e N and at least one n-i > 1 (see e.g. [6,7,8,9,10,11,12,13] and refs. 
therein). They are characterized by their genus number g, which is equal to 1 
for integrable systems and oo for chaotic systems. For pseudointegrable sys- 
tems, 1 < g < oo. In the specific case of systems with only right angles, 
g = 1 + iV 37r / 2 , where N 3n / 2 is the number of angles 3n/2. The different univer- 
sality classes persist when the systems are considered as quantum mechanical 
potential wells. In this case, the statistics of the eigenvalues shows characteris- 
tic features that correspond to the different dynamics of the classical systems. 



2 Spectral Statistics 

The first way to calculate the spectral statistics starts with the distribution 
P(s) of the normalized distances = (s i+ i — £i)/{s) between two consecutive 
eigenvalues £j+i and with the mean distance (s). For integrable systems, 
P(s) follows the Poisson distribution, whereas the of chaotic systems are 
Wigner-distributed [14,15]. Here, we consider the spectral rigidity A 3 (L) [16], 
which describes the mean correlations in a normalized energy interval of length 
L. We start from the integrated density of states N(e) of the normalized 
("unfolded") e where (s) = 1, a staircase function that can be approximated 
by a straight line. A 3 (L) is defined as the least square deviation between N(e) 
and its best linear fit r\ — r 2 e 



in average over many midpoints Eq from an interval Eq G AE = [0.5, 1.5]. The 
constants have been set to H — 2m — a — 1, where a is the lattice constant 
and m the mass, leading to an upper band edge -E max = 8. The range of Eq 
corresponds to wavelengths being 5 — 9a, much smaller than the size of the 
smallest boundary "teeth" that scatter around 16a for the v = 2-drum and 
60a for the v = 1-drum. AE contains 3000 — 5000 values of e, whereas the 
integration range L contains only up to 500 values (see below). The limit- 
ing curves for not too large L are A 3 (L) = L/15 for integrable systems and 
A 3 (L) = ln(L)/7r 2 — 0.07/7r 2 + O(L _1 ) for the ensemble of Gaussian orthogonal 
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matrices (GOE) [14,15], which serves as a good generally accepted limit for 
chaotic systems. For large L, A 3 reaches a plateau that depends on the small 
orbits. 



The spectral rigidity A 3 (L) can also be calculated by periodic orbit theory 
(see Fig. 1 for some examples of periodic orbits). In pseudointegrable systems, 
all periodic orbits form families of equal lengths and stay stable when their 
starting points are shifted to at least one direction along the boundary. The 
simplest families are the "neutral orbits" that bounce between two parallel 
walls. Using semiclassical quantum mechanics and neglecting additional con- 
tributions coming e.g. from diffractive orbits (that start and end at a salient 
corner), A 3 (L) under Neumann boundary conditions is given by [10,11,15], 
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The double sum is carried out over all orbit families (including repetitions, 
but only in forward direction) of lengths £i and £j. di and dj are the areas 
in phase space that are occupied by the respective orbit families and Hij = 
F(yi-yj) - F(yi)F(yj) -3(dF(y i )/dy i )(dF(y j )/dy j ), where F{y) = (siny)/y. 
The argument L enters via yi = (LiiA) / (IGiry/E^) with the system area A. 
It has been shown in [13] that for pseudointegrable systems, the double sum 
can be restricted to the diagonal terms of % = j and very long orbits can be 
neglected. We therefore need the and Oj up to some maximum length. 



3 Results and Conclusions 



We consider the systems shown in Fig. 1, i.e. two fractal drums of v — 1 and 
2, and the simple two-step system that is shown for comparison. First, we de- 
termined A 3 (L) by Eq. (1) from the eigenvalues calculated numerically by the 
Lanczos algorithm. Second, we calculated the periodic orbits by the hypothet- 
ical orbit method [12,13] that uses the fact that the systems are constituted 
by segments of lengths Xj and Yj in x- and y-direction, respectively (compare 
Fig. 1(b)). All periodic orbits have to pass through the segments an integer 
number of times. Therefore, only certain discrete angles ip (between trajectory 
and boundary) and lengths £ v are possible, where tany? = Yl,i n i^il J2j m jXj, 

?tp — 2 (Y.i'iT'iYi) 2 + (YsjmjXj) 2 , with positive integers and rrij. This 
restricts the number of possible orbits to certain pairs (<p(rii, m*), £(rii, m^)), 
which have to be tested. Due to the shielding role of the corners, not all 
hypothetical orbits really occur in each system. Therefore, the hypothetical 
orbit method checks which trajectories actually return to their starting point 
within the correct length of the trajectory. Unfortunately, the number of 
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Fig. 1. Geometry and periodic orbits: (a) The fractal generator, (b) The L-shaped 
system of g = 2. (c) Fractal drum of v = 2 with the generator applied to the upper 
border and g = 28. (d) Fractal drum of v = 1 with the generator applied to the 
upper border and g = 4. The segments of the drums are slightly deformed in order 
to avoid degeneracies of the orbit lengths. Some selected periodic orbits are shown 
for each geometry and the beam-splitting property of the salient corners and the 
different segments X\, X2, Y\ and Y2 are demonstrated in (b). 

hypothetical orbits increases with the number of segments by a power-law and 
the method becomes very time-consuming for systems with many corners. 

As a test for the orbit families we use a sum rule for the number N(£) of orbits 
with lengths smaller than £ (proliferation rate) [10,11], 

m = ^L, (3) 



where (a(£)) = J2i,£i<£ a %l J2i,e t <e 1 is the average area in phase space, occupied 
by the orbits with lengths smaller £ and b is a constant, depending slightly 
on the details of the system. On systems with g < 3, it has been found 
that bo = 1/4 for integrable systems and slightly larger for pseudointegrable 
systems, whereas (a(£)) is saturating to (a(£)) ~ 4A in integrable systems and 
to a much smaller value for pseudointegrable systems [11,13]. In Fig. 2(a), 
we plot the normalized average area (a(£))/(4:A) versus £ for the systems of 
Fig. l(b-d) with g — 2, 4 and 28. The decrease of (a(£)) with g is very drastic, 
which shows clearly the beam-splitting property of the salient corners. Each 
orbit that is disturbed by a salient corner splits into two different families, 
each of them covering a smaller area than without this corner. Therefore, 
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Fig. 2. (a) The normalized average areas {a)/{AA) of the periodic orbit families, 
with the area A of the systems is plotted versus the orbit length £ for the systems 
of Fig. 1 with g = 2 (triangles), g = 4 (circles) and g = 28 (squares), (b) S(£) for 
the same systems is plotted versus £. 
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Fig. 3. (a) The proliferation rate N{£) is plotted versus £. The lines show the the- 
oretical curves, whereas the symbols indicate the number of orbits that have been 
found by our calculations, (b) A^(L) is plotted versus L. Filled symbols refer to pe- 
riodic-orbit calculations (Eq. (2)) and open symbols to calculations from the eigen- 
values (Eq. (1)). (Same symbols as in Fig. 2.) The expected curves for integrable 
and chaotic systems are indicated by a solid and a dashed line, respectively. 

(a{£)) decreases with g, whereas N(£) increases. In Fig. 2(b), we plot S{£) for 
the same systems and determine the constant b by the slopes. We can see 
that b is very close to 0.25 (the value of integrable systems) for the system 
with g = 2 and 4 and considerably larger for g = 28. 

With the values of &o and (a), we can test the proliferation law N(£) of Eq. (3) 
for each system (see Fig. 3(a)). The lines show the theoretical curves, whereas 
the symbols indicate the number N of orbits found by our calculations. They 
agree very well for g = 2 and 4, whereas deviations from the quadratic shape 
of N(£) occur for g = 28, showing that many comparably small orbits are 
still missing in the spectrum of the v = 2-system. Due to the three different 
length scales of this drum and to the large number of orbit families, it is very 
time-consuming to find a sufficient number of them. However, we can already 
see by the present data that N(£) increases drastically with g. 
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We finally calculate A^(L) (i) from Eq. (1) using the eigenvalues and (ii) - 
for the systems of g = 2 and 4, where N(£) has been reasonably well approxi- 
mated by the calculated orbits - by periodic orbit theory. In both cases, the 
data from the eigenvalue and from the periodic orbit calculations agree very 
well. This indicates that higher-order terms, arising e.g. from diffractive orbits 
(additional to the periodic orbits) seem to be small (at most in the order of 
magnitude of the error bars arising from the Lanczos algorithm and the dis- 
cretization of the system). It will be very interesting to see, if this stays true for 
the fractal drum of v — 2, where many more corners and therefore diffractive 
terms are present. We also see in Fig. 3(b) that the curves decrease rapidly 
with increasing g, thereby coming closer to the Wigner curve for chaotic sys- 
tems. It will be very interesting to see, if drums of higher fractal generations 
v can come arbitrarily close to chaotic systems and in which way lengths and 
areas of the periodic orbits change with g — > oo. 
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